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Abstract 

Despite its appeal, real and simulated glass forming systems do not undergo an ergodic- 
nonergodic (ENE) transition. We reconsider whether the fluctuating nonlinear hydrodynamics 
(FNH) model for this system, introduced by us in 1986, supports an ENE transition. Using non- 
perturbative arguments, with no reference to the hydrodynamic regime, we show that the FNH 
model does not support an ENE transition. Our results support the findings in the original paper. 
Assertions in the literature questioning the validity of the original work are shown to be in error. 







I. INTRODUCTION 



It is appealing to associate the vitrification of tlie liquid into a frozen glassy state as a 
ergodic nonergodic (ENE) transition. Unfortunately there is strong evidence against the 
ENE transition scenario in physical and numerical experiments. This is in agreement with 
the results we found twenty years ago in Ref. |l| ( hereafter mentioned as DM) when we 
introduced the model of fluctuating nonlinear hydrodynamics (FNH). We present here a 
nonperturbative analysis of the FNH model and the possibility of an ENE transition. In the 
end our results here agree with those in Ref. There is no sharp ENE transition in the 
FNH model. Recent reservations 2| concerning our results in are shown to be unfounded. 
We also address some misrepresentations!^ of our work. 

In the theory of Classical Liquids, a new approach to studying the complex behavior 
of the supercooled state started with the introduction of the self consistent mode-coupling 
theory (MCT) 0, . The model referred to here is based on a nonlinear feedback mechanism 
due to the coupling of the slowly decaying density fluctuations in the supercooled liquid. 
The feed back effects at metastable densities strongly enhance the transport properties of 
the liquid. In the simple version proposed initially [a, Q, H] a sharp ergodic to non-ergodic 
(ENE) transition of the liquid into a glassy phase was predicted. This transition occurs 
at a critical density ( or at the corresponding values of other controlling thermodynamic 
parameters) beyond which the density auto correlation function freezes at a nonzero value 
over long times. Soon afterward it was demonstrated that this sharp ENE transition is |l| 
rounded. The absence of a sharp ENE transition in the supercooled liquids was supported 
by work using similar theoretical models. Two recent works has called these 

conclusions into question. The purpose of the present paper is to show that our previous 
analysis withstands careful scrutiny and to reassert that the results of ref. |l| are correct 
and captures the right phenomena for the removal of the ENE transition. 

We organize this paper as follows. In the next section we briefly introduce the FNH 
model. This is followed by an analysis of whether this model supports an ENE transition. 
In section HI we compare our findings here to those in DM. Next we comment on the works 
which question the conclusions in DM. We end the paper with a short discussion. 
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II. THE FLUCTUATING NONLINEAR HYDRODYANMIC MODEL 

In Ref. [l| a model for the long time relaxation behavior of the supercooled liquid was 
constructed using fluctuating nonlinear hydrodynamics. The dynamics of collective modes 
in the liquid was formulated with nonlinear Langevin equations involving bare transport 
coefficients. These nonlinear stochastic equations for the time evolution of the conserved 
densities are plausible generalizations of the macroscopic hydrodynamic laws. The set of 
collective variables {ipi} for the liquid we considered consists of mass and momentum densi- 



ties {p(r, t), g(r, t)}. The construction of the equations of motion ll| for the slow variables 
involve a driving free energy functional F which is expressed in terms of the hydrodynamic 
fields, i.e., p and g. The corresponding equilibrium distribution for the system is exp{—/3F). 
The free energy functional F is separated in two parts, F = F k\s , p] + Fu[p]. The depen- 



dence of F on g is entirely in the kinetic part Fx in the form 12| constrained by galilean 
invariance: 

The potential part F[/ is treated as a functional of the density only. The density p follows 
the continuity equation 

| + V,g^O. (2) 

having the fiux as the momentum density g which itself is a conserved property. The 
nonlinear equation for the momentum current density Qi is a generalized form of the Navier- 
Stokes equation 

f = .^V,lf]-pv/-^-ZLtA^O,. (3) 

The noise 6i is assumed to be Gaussian following the fiuctuation dissipation relation to the 
bare damping matrix L°-. For compressible liquids, the 1/p non-linearity appear in two 
terms in the generalized Navier-Stokes equation. These are respectively the convective term 
coupling two fiow fields and the dissipative term involving the bare viscosity of the liquid. 
The appearance of this non linearity in the hydrodynamic equations is formally avoided in 
Ref. [ij by introducing the local velocity field V(x, t), 

g(x,t) = p(x,t)V(x,t). (4) 
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The set of fluctuating variables in terms of which the renormalized field theory is constructed 
in our analysis therefore consists of the set ipi = {p, g, V}. 

The consequences of the nonlinearities in the equations of motion, i.e., renormalization 



of bare transport coefficients, are obtained using graphical methods of field theory [13| . The 
correlation of the hydrodynamic fields involve averages defined in terms of the action A 
which is a functional of the field variables {ipi} and the corresponding conjugate hatted 
fields {ipi} introduced in the MSR formalism. Using the equations of motions ([2]) and ([3]) 
respectively for p and g the action functional is obtained as [ll. 



. + ip 



9p „ 



(5) 



The theory is developed in terms of the correlation functions. 



G„^(12) = (^^(2)^„(1)) 



(6) 



and the response functions, 

G„^(12) = (^^(2)V^,(1)) . (7) 

The averages here are functional integrals over all the fields weighted by e~^. 
The nonlinearities in the equations of motion ([3]) and (jl]) give rise non-gaussian terms in the 
action (JSj) involving products of three or more field variables. The role of the non gaussian 
parts of the action A on the correlation functions are quantified in terms of the self energy 
matrix which show up in the equation satisfied by the response functions and that satisfied 
by the correlation functions. We begin with the response functions which satisfy: 



{G^%,{13) - S^,(13)l G .(32) = 5(12)5 



(8) 



with self energies S^^ which can be expressed in perturbation theory in terms of the two- 
point correlation and response functions. Using the explicit polynomial form of the action 
([5]), the response functions are expressed in the general form, 



G 



N, 



D 



(9) 
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where the matrix is given in table I and the determinant D in the denominator is given 
by 



D = Pl(c^' - g'c^) + iL{u + iq^-i) . (10) 

The various quantities are defined such that pl and L are identified as the corresponding 
renormahzed quantities respectively for the bare density po, speed of sound squared Cg and 
longitudinal viscosity Lq. We have in terms of single-hatted or response self-energies:, 



PL = pQ-iT.yy (11) 

L = Lo + iJ^gv (12) 
gc^ = qcl + Egp (13) 

and 7 is defined in terms of the self energy element = q'j. One can also show that the 
correlation functions of the physical un-hatted field variables are given by, 

Gap = — ^ Ga^iiCf^uGyp (14) 

where Greek letter subscripts take values p, g, V, and the self energy matrix Cf^o is given by, 

Cfiu = 2/5 ^ Lq5 fii,^ fj_,g — ^fio ■ (15) 

The double-hatted self-energies S^i> vanish if either index corresponds to the density. This 
model does not have a complete set of FDR linearly relating correlation and response func- 
tions. 

However, using the time translational invariance properties of the action ([5]), we ob- 
tained in DM the following fiuctuation dissipation relation between correlation and response 
functions involving the field g in the form : 

GvAq, ^) = -2p-'lmGg^a{q, u) (16) 
where a indicates any of the fields {p, g, V}. 
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III. ERGODIC-NONERGODIC TRANSITION AND FNH 



Does this model have an ENE transition? To answer this question we first pose the 
conditions for such a transition. Suppose, due to a nonlinear feedback mechanism, the 
self-energy blows up at small frequencies: 

i:^g = -A5{uj) . (17) 

This is presumed to result from a persistent time dependence of the density correlation 
function. This hypothesis is motivated by the one-loop contribution and the physics of the 
viscosity blowing up as one enters the glass. Is this assumption compatible with the set 
of Dyson equations? What we mean by a nonergodic phase is that Gpp shows a 5-function 
peak at zero frequency. Putting Eq. (fT7|) back into Eq. (fT^ we obtain a 5{uj) peak in Gpp 
as long as the response function Gp^ is not zero in the u; — > limit. We assume, with no 
reason to expect otherwise, that the — limits of pl, 7, and L are nonzero. With these 
assumptions D is not infinite in the low frequency limit and Gpg^ and G\rg are nonzero in 
the low frequency limit. Then from Eq. ljTll) we find that Gpp, Gpv, and Gyv show a 6{u!) 
component. Since Ggg vanishes as — > as long as D{uj = 0) 7^ 0, the correlation functions 
involving a momentum density index do not show a 5-function peak at zero frequency. So 
it is necessary for an ENE transition that Gpg not vanish as c<j — 0. This requires that pl 
goes to a nonzero value in the zero frequency limit and the determinant D not blow up as 

If, as expected, the self-energy contibution 'y{uj = 0) 7^ then the correlation functions 
Gpv, and Gyv show a 6{lj) component. Now we apply the FDT (|T6l) . Since Gyp and Gyy 
blow up, it then follows from the FDT that the imaginary parts of the response functions 
Ggp and Ggy also blow up. However we also require simultaneously that D*D is bounded, 
and imaginary parts of both pi^qD* and {uo + iq^'-))D* diverge. But since both D' and D" 
denoting the real and imaginary parts of D are bounded so pl and 7 must diverge. However, 
if these quantities blow up then from ( JTOl) it follows that D must also blow up and we have a 
contradiction. The obvious conclusion is that the original assumption of a nonergodic phase 
is not supported in the model. The key self-energy contribution is 7. If for some reason 
this quantity vanishes at zero frequency then G^y, and Gyy vanish as u goes to zero. Then 
Gpy, and Gyy do not show a 5{ijj) component and one does not have the constraints on pl. 
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7, and D. In this case one may have an ENE transition in this model. 



IV. RELATION TO DM RESULTS 

The argument we give in the hydrodynamic regime in Ref. [l| is completely consistent 
with the results presented above. The simplest way of understanding the argument in the 
previous section is to look at the response function 

r = ^/^L + tL 

The renormalization of the longitudinal viscosity L is computed, see Eq. (1121) in terms of the 
longitudinal part T,^y of the corresponding self-energy matrix Sg^v, of the isotropic liquid, 



L{q,z)=Lo + ^^y{q,z). (19) 

If we ignore the self energy Sy^, the expression f|T8|) is identical to the conventional ex- 
pression for the density correlation function with the generalized memory function or the 
renormalized transport coefficient L{q,z). The dependence of Gpp on the self energy 
in the renormalized theory is a consequence of the non-linear term involving the V field in 
the MSR action ([5]) and is originating from the nonlinear constraint (jl]) introduced to deal 
with the 1/p nonlinearity in the hydrodynamic equations. Analyzing the expression f|T^ for 
the correlation functions and the FDT relation (fT6ll we obtain in the hydrodynamic limit 
the following nonperturbative relation between the two types of self energies contributing 
alternatively to the renormalization of the longitudinal viscosity, 

7^^(0, 0) = 2/3-1 [7^^ (0, 0) + lim (7^(0, a;)/.;)] (20) 

where we have used in the above following definitions, in the isotropic limit, S^- ~ —q'^'jgg, 
~ —il^lgVi and Sp^ ~ 11 pg- The relation (120|) which is obtained from the FDT 
relation f|T6l) only, implies that both the self energies and S^y have the same diverging 
contribution in the low frequency limit. In the simplified model it is this contribution in 
terms of density correlation function which constitute the feed back mechanism of MCT 
and leads to the dynamic transition beyond a critical density. The singular contribution to 
the renormalized transport coefficient L in (fT8|) is now obtained in terms of the self energy 
S^y. As a consequence of ( 120|) it also follows that the response function Gpp is equal to the 



corresponding density correlation function Gpp in the hydrodynamic limit. It is important to 



note here that ( contrary to the assertion in Ref. [2|) this relation is not forced by us, rather 
it follows as a natural consequence of the relations (1201) linking the response to correlation 
self energies. 

The asymptotic behavior of the density correlation function is inferred from Gpp. The 
denominator of (fTSi) for the response functions contain the self energy which in this case 
is crucial for the long time dynamics and understanding how the ENE transition is cutoff. 
The density correlation function ( in the small g, uj limit) only freeze due to the feed back 
mechanism i/the self energy matrix elements is zero - a result obtained in the earlier 
section. In this regard it is useful to note that for the — limit the quantity L{uj + 27(7^) 
in D does not diverge even when L ~ l/cj is getting large, since L'^q^ remains finite in the 
non hydrodynamic regime uj ~ q^. To leading order in wave numbers 0) = g^7, is 

expressed in terms of the self energy using the nonperturbative relation 

7^(0,0) = ^^7>(0,0), (21) 



where c is the sound speed introduced in (11311 . Note that the relation ( ]2T1) is also obtained 
from the same fluctuation-dissipation relation [W^] . 



V. ABL AND CR 

We now address the criticisms made in Ref. {2] on our work. ABL imply that we 
misapplied the FDT relating Gpp and Gpp in the hydrodynamical limit. These authors offer 
that we assumed a linear FDT from the beginning. DM clearly discusses the consequences 
of not having a complete set of FD relations. On the other hand the introduction of the 
9 = 6F/6p field by ABL for obtaining a FDT in the linear form has not yet been shown 
to be useful. The nonlinear contribution in 6 comes from the part 6Fk/Sp, which actually 
gives rise [ll to the term 'Vj{gigj/p) in the generalized Navier-Stokes equation. The latter 
is essentially the 1/p nonlinearity which we address in our model through the introduction 
of the variable V. In this regard we believe that the importance of linear FDT in the MSR 
formulation has been overemphasized by ABL. Indeed in the absense of a linear FDT the 
response functions lose their physical meaning and become mere computational tools. From 
a physical point of view however what is important is that the correlation functions are time 
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invariant w 
In Ref . 



lich is maintained as can be directly seen from the above equation ([T 



3|] the cut off mechanism of Ref. jj has been questioned by treating the highly 



nonlinear model described above using a rather naive approach. CR basically make some 



phenomenological manipulations on a Newtonian dynamics model [IJ] , ending with a mem- 
ory function description they claim, without proof, is related to our model. All subsequent 
discussion of our work made by these authors are based on this claim. Our model, as shown 
on examination of table I, satisfies at all stages the density conservation law. The memory 
function proposed in CR to represent our work, Eq.(5) there, does not satisfy this conser- 
vation law. Therefore the analysis of CR does not apply to the model we studied. None 
of their conclusions concerning our work have any validity. CR concede that there is no 
error in our calculation, rather they offer vaguely that our model itself is the problem! Our 
model, as shown on examination of table I, satisfies at all stages the density conservation 
law. The source of their error appears to be the naive assumption that this model can be 
represented in terms of a single memory function This work represents a fundamental 
misunderstanding of the problem. 

Though the authors of both papers, ABL and CR, seem to agree that finally the ENE 
transition does not survive they disagree with our analysis of the problem. The arguments 
put forward in Ref. [3| to rediscover that the transition is finally cutoff are rather vague 
and of descriptive nature. These authors only seem to conjecture that the transition will be 



cutoff nonperturbatively citing other recent works IGj 



VI. DISCUSSION 

The basic feedback mechanism of MCT is a consequence of simple quadratic nonlinearities 
in density fluctuations ( arising from purely dynamic origin) that is present in the pressure 
term of the generalized Navier-Stokes equation. The ergodicity restoring mechanism goes 
beyond this. The description in terms of coupling to currents is a physically appealing way 
of explaining the nature of the FNH equations ( expressed in a form which can be sensibly 
related to the hydrodynamics of liquids). It is in fact the full implications of the density 
nonlinearities in the dynamics that cuts off the sharp transition to nonergodicity. This is 



also reflected in the fact that the basic conclusions of Ref. pj] follow even if the re 



nonlinearity is considered in a different manner. In fact by formulating the model 171] only 



evant 
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in terms of the fields {p, g} tlie same conclusions implying the absence of the dynamic 



transition is reached as in Ref. 



The 1/p nonlinearity mentioned above is treated here in 
terms of a series of density nonlinearites. The self energy matrix elements Hgv and are 
absent from the theoretical formulation in this case and the cutoff kernel is obtained here 
from a different self energy element Hgp. 

Twenty years ago we had predicted that the feedback effects from mode-coupling of 
density fluctuations, when properly analyzed keeping consistency with concepts of basic 
hydrodynamics, results in a qualitative crossover in the dynamics. We presented here a 
selfcontained nonperturbative proof that FNH does not support an ENE transition. This 
new analysis is completely compatible with the results of DM, simulations and experiment. 
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